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Abstract 



1 Introduction 



Grand Unified Theories (GUTs) based on the gauge group 5*0(10) [HIS] have been investigated 
extensively. This choice of gauge group has several appealing features. First of all, it has 
room for a right-handed neutrino per generation in the 16-dimensional irreducible spinor 
representation which includes all known matter fields. Thus it provides a beautiful explanation 
of the smallness of the neutrino mass via the "seesaw mechanism" |3j . Moreover, the existence 
of very massive right-handed neutrinos might also allow to explain the asymmetry between 
matter and antimatter in the Universe by thermal leptogenesis Furthermore, 50(10) 
contains the "Pati-Salam" [S] group SU{4:)c x SU{2)l x SU{2)ji as subgroup, meaning that 
parity is preserved at high energy and broken spontaneously. 

On the other hand, the fact that the rank of 5*0(10) is five causes some complications. 
Recall that the rank of the Standard Model (SM) gauge group Gsm = 5'f/(3)cx5'[/(2)iX[/(l)y 
is only four. There are several ways of breaking 5*0(10) down to Gsm, depending on which 
representations of Higgs fields are introduced in the theory. Here we consider the possibility 
of having intermediate phase(s) at energy scales well below the GUT scale. The existence of a 
scale near 10^^ GeV can be motivated by neutrino oscillation experiments O [TJ E]: the mass 
of the heaviest neutrino cannot be less than ^/^w^^ ~ 0.04 eV. In the seesaw mechanism 
this translates into an upper bound on the right-handed Majorana neutrinos mass if we 
assume that the largest neutrino Yukawa coupling is order unity. Mat < 10^"^ GeV. Note that 
Mn breaks the SU{2)r subgroup of 5*0(10). It thus seems natural to assume the left-right 
symmetric subgroup of 5*0(10) to be broken to Gsm near this scale ("Mr"), if we assume that 
the Yukawa coupling that gives rise to the Majorana mass M^v is also of order unity. 

In this work, we will analyze the consequences of this assumption, by considering the low 
energy phenomenology of a supersymmetric 5*0(10) model suggested by Aulakh et al. [9]. It 
features the symmetry breaking chain 

50(10) ^ G422D ^ O3122 OsM . (1) 

Mx Mc Mj{ 

Here we have used the notation O3122 = SU{3)c x U{1)b-l x SU{2)ji x SU{2)l and G422D = 
SU{4)c X SU{2)r X SU{2)l x D, where D is a discrete symmetry which ensures that SU{2)l 
and SU{2)r have equal gauge couphngs. We assume universal ("mSUGRA" [lOj) boundary 
conditions for the soft supersymmetry breaking terms at the GUT scale Mx- This means that 
all soft breaking scalar masses are equal to mo at the GUT scale, while all gaugino masses are 
equal to M1/2; moreover, all SUSY breaking trilinear scalar couplings are characterized by the 
single parameter Aq. 

Introducing two intermediate scales, and the corresponding additional gauge, matter and 
Higgs superfields, has three main effects. First, the right-handed neutrinos obtain Majorana 
masses at scale M/j by coupling to the 126-dimensional Higgs whose vacuum expectation 
value (VEV) is responsible for breaking G3122 in Eq-dlD- These Majorana Yukawa couplings, 
as well as the extra Dirac couplings of the light neutrinos, will change the low energy spectrum 
of soft breaking parameters via renormalization group equations (RGBs). Secondly, since we 
have to introduce many more additional Higgs than gauge superfields to achieve the symmetry 
breaking chain ([1]), all gauge couplings increase quite rapidly at high energy scales > Mr. As 
a result the gaugino masses, which we assume to be universal at Mx, decrease significantly 
when they evolve down to M/j. Finally, the enhanced gauge symmetry at energies > M^i also 
increases the size of gauge contributions to the RGE of all scalar masses. Note that the second 



and third effect tend to cancel, if tlie scalar masses are expressed in terms of the GUT-scale 
input parameter ttlq and Mi/2- 

The rest of this paper is organized as follows. In the next Section, we review the main 
features of the model jilj we are considering. We also describe the numerical methods used 
in our analysis. In Sec. 3 we discuss the most important experimental and cosmological 
constraints on the parameter space of the model. Our numerical results are given in Sec. 4. 
Special attention is devoted to the regions of parameter space where the lightest neutralino 
makes a good thermal Dark Matter candidate in standard cosmology. Finally, we conclude in 
Sec. 5. 



2 The Set-Up 
2.1 The model 

We will consider the model suggested by Aulakh et al. [9]. It is based on the gauge group 
5*0(10). Besides three generations of matter superfields residing in 16-dimensional represen- 
tations as well as the 45-dimensional gauge superfields, we introduce Higgs superfields in the 
54, 45, 126, 126 and 10 representations of SO (10). The Higgs superfields required to break 
SO (10) down to GsM can be described by the tensors 

54 Sij = ^ji ^-^id S^i = , 45 i A-ij —— 7 

126 . ^ijififji '^(^ijklmopqrs^ opqrs i 

i 



126 . Sjj/j;^ ^^(^ijklmopqrs^ opqrs (2) 

where the subscripts i,j, k, . . . run from 1 to 10, and repeated subscripts are summed. 

This allows us to realize the symmetry breaking chain ([1]) with a purely renormalizable 
superpotential, given by |9j 

WssB = ^tr^2 ^ ^^^^3 ^ !^trA2 + XtiA'S + m^SS + r]sJ:^S + rj^s's + r^^SSA . (3) 

^ O ^ 

A crucial observation [HI |9] is that some components of the Higgs superfields listed in (j2]) 
are much lighter than one might naively expect. For example, even though the 45-plet A 
is responsible for the breaking of G422D to (^3122 at scale Mc, some components of A only 
acquire masses of order M^/Mx or M'^/Mq, whichever is larger. Similarly, even though S 
and S are responsible for breaking G3122 to the SM gauge group, some of their components 
only get masses of order M'^/Mx- On the other hand, some components of A, S and E obtain 
masses of order Mx- 

This is summarized in Table [H Here we have used the decompositions of the Higgs fields 
under SU{A)c x SU{2)l x SU{2)r: 

S = (1,1, 1)0(20, 1,1) ©(1,3, 3) ©(6, 2, 2); 

A = (15, 1,1) ©(1,1, 3) ©(1,3,1) ©(6, 2, 2); 

S = (10, 1,3) ©(To, 3,1) ©(15, 2, 2) ©(6, 1,1); 

E = (To, 1,3) ©(10, 3,1) ©(15, 2, 2) ©(6, 1,1). (4) 



The components of the Higgs fields that acquire large vacuum expectation values (vevs) appear 
as the first term in each right-hand side (rhs) of Eqs.(jl]); in addition, the (1, 1, 3) component 
of A is also assumed to obtain a nonzero vev [9J. 





iVlcLoo 


all of S 

all of A, except (15, 1, 

all of S and S, except SU{4)c (anti-)decuplets 


~ Mx 


(10,3,1)2 and (10,3,l)s 

color triplets and sextets of (10, 1, 3)^ and (10, 1, 3)^ 
color triplets of (15, 1, 1)a 


~ Mc 


(5° -5'), 5+, r 


~ Mr 


color octet and singlet of (15, 1, 1)a 


-Mi^max 


(50 + 5'), 5++, 


~ M2 = Ml/Mx 



Table 1: The spectrum of Higgs superfields after symmetry breaking. The Higgs superfields 
have been introduced in Eq.([2]), and their decomposition into irreducible representations of 
SU{4)c X SU{2)l X SU{2)r is given in Eq.(jl]). form the color singlet part of the 

(10, 1, 3) component of S, while 5°'^' form the color singlet part of (10, 1, 3) of S. Adapted 
from ref . [9] . 

We also need Higgs superfields in the 10-dimensional representation of 5*0(10) to provide 
the Higgs doublet superfields of the Minimal Supersymmetric Standard Model (MSSM) that 
break the electroweak gauge symmetry. Minimal 5*0(10), with a single 10, would require 
all Yukawa couplings of one generation to unify, which leads to wrong predictions for ratios 
of quark and lepton masseso Introducing G422D as symmetry group between Mx and Mc 
aggravates this problem, since it predicts Yukawa unification at scale Mc if both MSSM Higgs 
doublets reside in a single (1,2,2) of SU{A)c x SU{2)l x SU{2)ji. We therefore include two 
such superfields. We assume that the additional bidoublet obtains a mass through the coupling 
to the (1, 3, 1) of A, in which case its mass will be of order M2 = M'j^/Mx [95- 

Let us discuss the structure of the matter Yukawa couplings in a bit more detail. Here 
we are only interested in third generation couplings, which can be large enough to affect the 
weak-scale sparticle spectrum significantly. The Yukawa unification conditions we will derive 
will not work for first and second generation fermions. We assume that this problem is solved 
by introducing some more complicated fiavor structures, e.g. via non-renormalizable terms, 
without introducing additional large couplings. 

At energy scales below M2 we have the well-known MSSM superpotential, 

WVuk,MSSM = YMVHu + YdD'QHd + YeE'LHd , (5) 

where Q and L are the quark and lepton doublets, V^, D'^ and E'^ the corresponding singlets, 
and Hu and the two Higgs doublet superfields. In Eq. ^ we have suppressed all generation 
and group indices. 



*This prediction can be made to work for the third generation, if the ratio of MSSM Higgs vevs tan (3 is 
large and sfermion masses he weh above a TeV [12]; however, they wiU fail for the first two generations. 



At energies above M2 the second pair of Higgs doublets as well as some parts of the SU{2) 
triplet Higgs superfields (see Table 1) become accessible. A general ansatz for the matter 
superpotential is then 

i=l 

The last term in Eq.([6]) results from the interaction giving rise to large Majorana masses 
for the right-handed neutrino superfields (see below). The light Higgs doublets H^, are 
mixtures of the Higgs superfields appearing in Eq.([n]): 

= cos iPuHu,i + sin v?«i^«,2 ; 
Hd = cos ipdHd,i + sin ipdHd,2 ■ (7) 

At scales above Mr, and D'^ form a doublet Q'^ of SU(2)fi.; similarly, the right-handed 
neutrino superfield iV^ and E'^ form an SU{2)ji doublet L^l^ Moreover, the Higgs superfields 
Hu,i, Hd^i are grouped into two bidoublets $j. Finally, at this scale all members of the SU{2)r 
triplet Higgs superfield 5 become accessible. The superpotential (jH]) then becomes 

2 1 

W^Yuk,3122 = (Y^^^Q'Q^^ + yl,^L''L^^) + l^Y^L'^L' . (8) 

i=l 

The last term in Eq.(l8]) gives rise to large Majorana masses for the A^'^ once the neutral 
component of the SU{2)r triplet 5 e S gets a vev. Finally, at scales above Mc, Q and L are 
unified into F in the (4,2,1) representation of G422, while and join to form F'^ in the 
(4, 1, 2) representation. One is then left with a single Yukawa coupling per Higgs bidoublet, 

2 1 _ _ 

W^Yuk,422 = ^^F'^F^, + -Yj, {F'T.rF' + FT.lF) , (9) 
1=1 

where T^r and are in the (10, 1, 3) and (10, 3, 1) representation, respectively, of S'[/(4)c x 
SU(2)l X SU(2)r; the last term in Eq.(l9]) also has to have coupling Yn due to the discrete D 
symmetry. 

As a first simplification, let us work in the basis where Y2 = 0. This can always be accom- 
plished by a unitary rotation between the two Since superpotential couplings renormalize 
multiplicatively, this choice is renormalization scale invariant. It is then easy to see that, 
through the matching conditions at Mc, Yq 2 = = in Eq.([H]); similarly, matching at scale 
Mr implies Yu^2 = ^d,2 = ^e,2 = in Eq.([2]). The sums in Eqs.dS]) and ([HD thus also collapse 
to single terms. Inserting Eqs.([7]) into Eq.([6]) then leads to the following matching conditions 
for the MSSM Yukawa couplings at scale M2: 

Yu,i = Yu/ cos ; Yd^i = Yd/ cos (pd ; ^e.i = YJ cos (pd ■ (10) 

We can get phenomenologically acceptable couplings only if (fu 7^ fd- 

Note that the high-scale couplings Yf^i are always larger than or equal to the low-scale 
(MSSM) couplings Yf (/ = u,d,e). On the other hand, we know that in the MSSM the top 



'''Note that Q'^ and L'^ are independent left-chiral superfields, not the charge conjugates of Q and L. 



Yukawa coupling is already fairly close to its upper bound imposed by the requirement that it 
remains perturbative up to very large scales. Eq. ffTOl) therefore implies that | cos ipu\ — 1- For 
definiteness we therefore set 

cosipu = l, (11) 

i.e. (fu = 0. This minimizes YuX, we will see shortly that it also minimizes all other MSSM 
matter Yukawa couplings above scale M2. These couplings appear with positive signs on the 
right-hand side of the RGE for the new coupling Y^- The choice ffTTj) therefore maximizes the 
upper bound on YnIMpi) that can be derived from the requirement that this coupling remains 
perturbative up to Mx- We will see below that this in turn minimizes the lower bound on the 
mass of the light physical neutrino for fixed M^. 

Eq.dHD implies that F«,i(M^) = 1^,1 (M^) = Yg^i{Mji). This is compatible with Eqs.([inD 
and (ITT]) only for 



(12) 



Yu{M2) 

the last factor in Eq. f|T2|) accounts for the different RGE running of Y^^i and Y^^i caused by 
the different hypercharges of the f/^ and Z)^ superfields. Since this factor is quite close to 
unity, Eqs. fllOI) - ffT^ imply that Y^^i and i are very similar. Since even in the MSSM the 
bottom and tau Yukawa couplings become similar at large energy scales, Eq. (fT2|) implies that 
all third generation Yukawa couplings will be comparable to the top Yukawa coupling at all 
scales above M2. In the given framework this is inescapable, unless we introduce additional 
heavy superfields which mix with the MSSM matter fields. 

At the SU{4:)c breaking scale Yg^i and Y/ i are unified into the single coupling Y^. The 
unification of the bottom and top coupling can always be achieved through an appropriate 
choice of (pd] however, the unification of the bottom and tau Yukawa couplings is a nontriv- 
ial constraint. This prediction is similar to that of minimal SU{5). In a scenario without 
intermediate scales, the tau Yukawa coupling at scale Mx is typically a bit larger than the 
bottom coupling. In our case unification should happen at Mc < Mx, which reduces the 
difference between the two couplings at their putative unification scale. On the other hand, 
above the scale Mi < Mc the SU{3)c coupling is larger in our scenario than in the MSSM. 
This increases the RG running of Yi,. The two effects largely cancel. As a result, we find that 
Yt-{Mc) exceeds Yi,{Mc) by typically 10 to 20%. We blame this on threshold effects - Table 1 
shows that quite a few new fields attain masses of order Mc - and/or on the additional physics 
required to reproduce masses and mixing angles of the lighter SM fermions. As a practical 
matter, we set 

yi(M,) = MM£)±M^. (13) 

The superpotential (IE]) generates neutrino masses through the celebrated ("type I") see- 
saw formula [3]jil 

„ '"^ (>»<g°»' (14) 

Here a E (10, 1,3) G S is the neutral component of the SU{2)fi triplet Higgs bosonJl Note 
that the neutrino Dirac Yukawa coupling is related to that of charged leptons by SU{2)r, 

■i-Note that there is no quartic scalar couphng which could lead to a "type II" seesaw contribution. 
^The field listed in Table 1 is the physical remnant of a after 6*3122 — > Gsm symmetry breaking. 



which in turn is related to the (top) quark Yukawa couphng by SU (4)c symmetry, as described 
above. We assume here that the Ya? are (mildly) hierarchical, so that only the third generation 
coupling is large enough to effect the weak-scale spectrum significantly via the RGE. 

For given = (a) and light neutrino mass m^, Eq. flT^ can then be used to determine 
the value of Y^r at the GUT scale. We vary mj, between 0.2 and 0.4 eV. Note that smaller 
values of lead to a larger coupling Y^. 

The occurrence of fields that are not part of the MSSM at mass scales well below Mr is 
crucial. As well known, in the MSSM all three gauge couplings (almost) meet at an energy 
scale near 2 ■ 10^^ GeV [13]. Without additional fields that are lighter than M/j it would not 
be possible to modify the running of the gauge couplings such that intermediate scales, and 
hence energy ranges where the symmetry group is larger than Gsm but smaller than the GUT 
group, can occur. We will analyze the running of the gauge couplings in more detail in Sec. 4. 

We will see in Sec. 4 that the lightest new particles, with mass M2 ~ M'^/Mx-, are still much 
too heavy to directly lead to visible effects at collider or rare decay experiments. Nevertheless 
their existence affects the renormalization group equations (RGE) describing the running of the 
masses of all superparticles and Higgs bosons. The one-loop RGE for Yukawa couplings and 
soft breaking parameters that hold for different ranges of energies are listed in the Appendix. 
In order to compare with the frequently studied [H] mSUGRA or cMSSM scenario, we assume 
universal boundary conditions, as already noted in the Introduction. 

2.2 The numerical calculation 

The RGE listed in the Appendix are too complicated to allow an analytical solution. Instead, 
we incorporated them into the code S0FTSUSY2 . [15] . This program computes the weak-scale 
MSSM spectrum by iteratively solving the RGE, starting from universal boundary conditions 
for the soft breaking parameters. An iterative treatment is necessary since many parameters 
are fixed at the weak scale, rather than the GUT scale. These include the three (MS)SM 
gauge couplings, the masses of SM matter fermion^, the mass of the Z boson, and the 
ratio tan/? of vevs of the two MSSM Higgs bosons. We use one-loop RGE throughout, but 
include important weak-scale threshold corrections; these are known to change the physical 
masses of third generation fermions significantly, in particular at large tan/3 [TH]. Note that 
the program implements radiative breaking of the electroweak gauge symmetry [17], again 
including important weak-scale threshold corrections. 

At the intermediate scales Mi, M2, Mr and Mc (some of) the RGE have to be changed. 
In the discussion of Yukawa couplings we described how to pick the appropriate high-scale 
couplings, given the low-energy couplings. This procedure is applicable when going from low 
to high energies. When going in the opposite direction, we use the same matching conditions, 
employing the values of cosy^d and the ratio of Yi^i{Mc) /Yq^i{Mc) determined from the previ- 
ous RG running from low to high energies to fix the values of low-scale Yukawa couplings. The 
matching of gauge couplings and soft breaking terms directly follows from the group structure, 
and will be discussed in Sec. 4.1 and in the Appendix, respectively 

The output of SOFTSUSY is passed on to the program micrOMEGAs 1.3.7 [18j, which com- 
putes the Dark Matter (DM) relic density as well as the BR(6 — > 57) and 5a^, the anomalous 
magnetic moment of the muon (see below). 

^Only the masses of third generation fermions are kept, since the Yukawa couphngs in the first and second 
generation are too smaU to significantly affect the evolution of the sparticle masses. 



3 Accelerator and Cosmological Constraints 



In this Section we describe the constraints we impose on the model. 

3.1 Electroweak symmetry breaking and tachyons 

As mentioned earher, electroweak symmetry breaking (EWSB) is incorporated into SOFTSUSY. 
Technically, it solves equations that allow to express /i^ and the bilinear Higgs soft mass 
parameter Bfi in terms of Mz and the ratio of vevs tan (3. However, these equations sometimes 
formally lead to fi"^ < 0, which indicates that EWSB is not possible for the given set of input 
parameters. For reasons that will become clear shortly, here we are mostly interested in 
solutions with large tan /?. In this case EWSB is possible iff the (properly threshold corrected) 
value of the squared soft breaking mass of the up-type Higgs boson at the weak scale is 
negative, m^^(MsusY) < 0. 

3.2 Constraints from collider searches 

As in mSUGRA, the most important constraints are those on the masses of the lightest Higgs 
boson and the lightest chargino. In combination, they imply that constraints on the masses 
of strongly interacting sparticles [I9] are automatically satisfied. 

We interpret the limit Mhc.^ > 114.4 GeV, which comes from searches for e~^e~ ZH^, 
as imposing a lower mass on the mass of the lighter CP-even Higgs boson of the MSSM, 

Mft>lllGeV, (15) 

where we allowed for a ~ 3 GeV theoretical uncertainty [20] in the calculation of m^. We also 
require 

m-±>104GeV, (16) 

since scenarios allowing chargino masses significantly below the highest LEP beam energy 
cannot be realized in our scenario: these scenarios all require the presence of sneutrinos with 
mass near or slightly below that of the chargino, and scenarios where both the sneutrino and 
xf have mass below the limit (fT6|) violate the Higgs constraint (fT5|) . 

3.3 Branching ratio of 6 ^ 57 

In the SM, flavor changing neutral currents (FCNC) are absent at tree level. Thus, the 
radiative B — > decay is mediated by loops containing up-type quarks and W bosons. As 
well known |21] , SUSY loop contributions can be comparable to those from the SM. Therefore, 
the measurement of the branching ratio for this decay, performed by CLEO, Belle and BaBar 

B{h 57) = (355 ± 24+?o ± 3) X 10"^ (17) 

can be used to constrain the parameter space of our model. The first error in f|T7j) includes 
statistical, systematic, extrapolation and b ^ d'j contamination errors, while the last two are 
estimated to be the difference of the average after varying the central value of each experimental 
result by ±la. To be conservative, we take the linear sum of the errors, since the calculation 
strongly depends on the assumptions of the boundary conditions. Even minor deviations from 



strict universality, for example due to the running between Mx and Mpi [23l [21], can have 
very large effects [25] . 

As mentioned above, we used micrOMEGAs 1.3 [T8] to calculate the branching ratio. 
Therein, minimal flavor violation (i.e. the only source of flavor violation at the weak scale is 
in the CKM matrix) is assumed ^26j; hence only contributions from charged Higgs and top 
quarks, and charginos and stops are included. These contributions are indeed usually by far 
the dominant ones if universal boundary conditions are assumed [27], as in our analysis. 

3.4 The anomalous magnetic moment of the muon 

The anomalous magnetic moment of the muon is one of the most precisely calculated and 
measured quantities. There is an about Sa discrepancy between the SM prediction based 
on data from e^e~ annihilation into hadrons and the experimental value. While this is still 
somewhat controversial - an SM prediction which instead makes use of r decay data plus some 
assumptions is in fair agreement with the data - we here want to investigate the parameter 
space of our model that allows to explain this discrepancy. 

The world average, dominated by data from the E821 collaboration at BNL, is [19] 

^exp ^ ^ (1165920.80 ± 0.63) x 10"^ . (18) 

The theoretical value [2S] is calculated as the sum of (i) pure QED contributions including the 
diagrams of virtual photon, vacuum polarization (VP) from e, fi and r, and leptonic light-by- 
light scattering, (ii) hadronic contributions including VP from quarks, most reliably estimated 
using e~^e~ — > hadrons data, and hadronic light-by-light scattering, and (iii) electroweak 
contributions. The resulting SM prediction is [19\ 

^theory ^ (1165919.52 ± 0.52) X 10"^ . (19) 

Demanding that supersymmetric loops, involving smuons and neutralinos or smuon neutrinos 
and charginos, lead to agreement between theory and experiment at the 2a level thus implies 

4.6 X 10"^° < 5a^,susY < 21.0 x 10"^° . (20) 

We use micrOMEGAs to calculate 5a^,susY- 

3.5 Dark Matter reUc density 

We assume that all cosmological Dark Matter (DM) consists of lightest neutralinos. This 
implies that Xi has to be the lightest superparticle (LSP); this imposes a constraint on the 
parameter space of our model. 

Far more important is the requirement that the thermal Xi relic density, calculated using 
micrOMEGAs under the usual assumptions of the minimal cosmological model [21], reproduces 
the value derived from the WMAP 3-year data [30j and other observations pertaining to 
structure formation in the universe: 

0.097 < n^Mh^ < 0.113, at 68% CL . (21) 

Here ^dm is the DM mass density in units of the critical (closure) density, and h is the Hubble 
constant in units of 100 km/(s-Mpc). As we will see below, this provides the most stringent 



constraint on the parameter space. This is not surprising, given the small size of the error 
bars in (|2T1)FI 



4 Results 

We are now ready to present some results. We begin with an analysis of the running of the 
gauge couplings, which determines the values of our intermediate scales. We then discuss 
analytical results for first and second generation sfermion as well as gaugino masses, before 
analyzing the ratios of (s)particle masses that are most relevant for the calculation of the DM 
relic density. We will conclude this Section with a survey of the parameter space of the model. 

For the top quark mass, we have taken = 170.9 GeV, as has recently been measured at 
the Tevatron [32] . 

4.1 RG Analysis of the gauge couplings 

The current world averages of the gauge coupling constants at scale Mz are 

ai{Mz) = 0.01695, a2{Mz) = 0.03382, a^iMz) = 0.1176. (22) 

Note that we use GUT normalization for the U{1) gauge coupling, i.e. our ai exceeds the 
hypercharge coupling ay in its usual normalization by a factor of 5/3. The values of these 
couplings at different energies are determined by RGE; to one-loop order, these can be written 
as 

^ = -f^&. (^ = 1,2,3). (23) 
at 2n 

Here t = ln[Q/Qo), where Qo is some reference energy scale. Note the minus sign in Eq. fl23l) : 
in this convention, a positive bi corresponds to an asymptotically free gauge coupling. The 
values of the bi depend on which particles are "active" at a given energy scale Q; in the usual 
step function approximation of integrating out heavy particles, we treat all particles with 
masses < Q to be (fully) active at scale Q. 

This leads to the values of the bi listed in Table [21 which we adapted from ref . [9] . Note that 
we list the coefficients that allow to describe the running of the three factor groups of the SM 
gauge group. The SU{2)l factor remains independent up to scale Mx, i.e. the third column 
of Table [2] always describes the running of the coupling of an SU{2) group. Recall that above 
M(7, SU{2)l and SU{2)r have the same coupling, since the discrete symmetry D is exact; 
the coefficient b'^^ therefore also describes the running of the SU{2)fi coupling. Moreover, 
at scale Mq the strong interactions get embedded into SU{4:)c, with boundary condition 
g3{Mc) = gi{Mc)- The coefficient b^^^ therefore describes the running of the SU{A)c gauge 
coupling, which is the same as the running of the coupling of the SU{?>)c subgroup of SU{4)c 
at Q > Mc. 

The fate of the U{1)y factor of Gsm is a bit more complicated. At scale Mr it gets 
embedded into SU {2)r x U{1)b~Li with matching condition a^^ = 2/ (5aB_L) + 3/ (5q;2_r). The 
U{1)b-l factor in turn gets absorbed into SU{4:)c at scale Mc, i.e. aB-L{Mc) = a^^Mc). 
Although the hypercharge coupling is thus "spread" over two different gauge couplings for 
Q > Mji, its running can still be described by Eq. (l23|l . with coefficient listed in Table [2l 

* Recently the WMAP 5 year data have been released [31]. The resulting range for ^l^uh^ is very similar 
to that of Eq. (PT|l . 



Table 2: The coefficients of the beta functions of the gauge couphngs of the SM gauge group, 
vahd at different energy scales Q. 



Energy range 
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Mz<Q<Ms 


-41/10 


19/6 
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Ms<Q <M2 


-33/5 


-1 
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M2<Q<Mr 


-12 


-2 
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Mr<Q < Ml 


-48/5 


-2 


3 


Mi<Q < Mc 


-48/5 


-2 





Mc<Q <Mx 


-194/5 


-42 


-34 



Eqs.( l22i) and ( l23l) . together with the coefficients h] listed in Table Ej allow us to predict 
the values of the gauge couplings at all Q > Mz- Of course, the three gauge couplings of 
the (MS)SM are supposed to meet at scale Mx in our model. This leads to two independent 
constraints. On the other hand, the intermediate scale Mr and Mq are free parameters of our 
model; the scales Mi and M2 are derived quantities, as described in Table [H For given value 
of Mx the two independent unification conditions can thus be solved for Mr and Mq- The 
running of any one of the three (MS)SM gauge couplings can then be used to determine the 
value of the 5'O(10) gauge coupling au- Notice that this procedure will work for any assumed 
value of Mx, i.e. it still leaves one parameter undetermined. We refer the reader to ref. [9] for 
a further discussion of the unification condition, including explicit solutions of the RGE of the 
gauge couplings. 

In Fig. [1] we show one-loop predictions for the intermediate scales Mr and Mc, as well as 
the value of Mx, as function of l/au- We see that smaller values of Mx correspond to larger 




Figure 1: The values of the intermediate scales Mr and Mc, as function of the inverse of the 
5*0(10) coupling au. Here we took Ms = 1 TeV as sparticle mass scale. 



values of au- The reason is that decreasing Mx increases the ratios Mx/Mr and Mx/Mq. 
Table [2] shows that all 6j are large and negative for Q > Mq] recall that this corresponds to 
gauge couplings increasing with energy. A large Mx/M^ means that these beta-functions are 
valid over a large range of energies, leading to a large value of ajj. On the other hand, proton 
decay through dimension 6 operators conservatively requires Mx > 3 ■ 10^^ GeV. Fig. [T] shows 
that this corresponds to au — 1/13.5, safely in the perturbative region (significantly smaller 
than az{Mz), for example). 

Since the purpose of our paper is to study the influence of the intermediate scales on 
the low-energy spectrum, we take this minimal value of Mx as our default choice. The 
intermediate scales are then found at 

Mr = 10^^-^^ GeV , Mc = lO^^'^^ Q^y _ (24) 

Increasing Mx reduces the impact of the intermediate scales. At Mx — 10^^'^ GeV, corre- 
sponding to au ~ 1/21, the scales Mr and Mq coincide. When Mx is increased to about 
10^^'^ GeV, Mc in turn coincides with Mx- At that point no intermediate scales are left, i.e. 
this limit reproduces the usual MSSM. Higher values of Mx are not possible. By varying Mx 
between 10^^'^ GeV and 10^^'® GeV we can thus smoothly turn on the intermediate scales and 
study their impact on weak-scale physics. 
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Figure 2: The running inverse gauge couplings a^'^. The curves at the top [middle, bottom] 
are for the f/(l)y [SU{2)l, SU{3)c] couplings. Solid and dashed curves show results for one- 
and two-loop RGE, respectively. 

The rapid increase of the gauge couplings at Q > Mq can also be seen in FigjJl which 
shows the running of the gauge couplings as function of the energy scale for our default set 
of parameters. The solid lines show the predictions from the one-loop RGE we have used 
so far, whereas the dashed curves are based on two-loop RGE [33] (ignoring, however, the 
subdominant contributions from Yukawa couplings to the running of the gauge couplings). 
Evidently using two-loop RGE increases the intermediate scales for this value of Mx, making 




the model more mSUGRA-like. However, an analysis based on two-loop RGE should also 
treat the (rather numerous, in our case) threshold corrections more carefully. So far we have 
assumed that all (s)particles whose masses are of the order of a given scale, as listed in 
Table [U have exactly that mass. This will not be true in many cases. However, the exact 
masses will depend on many unknown couplings describing interactions of these superheavy 
fields. A proper treatment of threshold corrections would therefore introduce many new free 
parameters. Since threshold and two-loop effects are generically of similar magnitude [33], we 
assume that there are combinations of parameters where an analysis including two-loop and 
threshold effects leads to similar results as the one-loop analysis. The use of one-loop beta 
functions has the practical advantage that the equations determining M/j and Mc can easily 
be solved analytically 

4.2 Analytical results 

At the one-loop level, the gaugino masses evolve in the same way as the squared gauge 
couplings do. Therefore, the ratios of weak-scale gaugino masses are the same as in mSUGRA, 
i.e. Ml : Ma : Ms ~ 1 : 2 : eEl This follows from the fact that the three MSSM gauge 
couplings are identical at Mx, and have their measured values (l22l) at scale Mz- These ratios 
are therefore independent of the intermediate scales. 

However, for fixed M1/2 the weak-scale gaugino masses are now much smaller than in 
mSUGRA, since the ratios ai{Mz)/au are much smaller, as shown in Fig. [21 Writing 

Mi(MsusY) = qMi/2 (z = 1 , 2 , 3) , (25) 

we have 

ci~0.23, C2~0.46, C3~1.4, (26) 

for MsusY ~ 1 TeV; these are nearly two times smaller than the corresponding coefficients in 
mSUGRA [lOj. 

The RGE for the masses of first and second generation sfermions, whose Yukawa couplings 
are negligible, can also be solved analytically [TD]. Writing 

^^(Msusy) =ml + cjMli^ , (27) 

we have 

Ce^ ~ 0.15 , - 0.21 , Cg - 1.16 . (28) 

Here cr and II stands for U{1)y singlet and doublet sleptons, respectively, while q stands for an 
average first or second generation squark; as in mSUGRA, SU (2) l doublet squarks are shghtly 
heavier than singlet squarks. We checked that the analytical and numerical calculations of 
match within 0.1%. Note that the coefficient cg^ is numerically almost the same as in 
mSUGRA [TD]. This is due to a cancellation of two effects. On the one hand, is a non- 
singlet under both SU(2)r and S't/(4)c giving rise to new gauge contributions to its mass at 
scales above Mr and Mc, respectively, which increase cg^. On the other hand, we saw that 

* These are running masses. The on-shell masses differ by weak-scale threshold corrections [IF, which are 
included in SOFTSUSY. 

^The running weak-scale sfermion masses also receive small Z?— term contributions, which we omit in the 
following discussion, but include in the numerical analysis. In addition, the physical (pole) masses again differ 
from the running masses by threshold corrections (16j, which are included in SOFTSUSY. 



for fixed M1/2 the gaugino masses at scales Q < Mx are smaller than in mSUGRA, which 
reduces all Cj. The latter effect is dominant for all fields that transform non-trivially under 
either SU{2)l or SU{?>)c- As a result, the mass difference between SU{2)l singlet and doublet 
sleptons is significantly smaller than in mSUGRA; recall that the SU{2)l doublet sleptons are 
singlets under SU{2)ji. 

Note that if we apply the universal boundary conditions at some energy scale Q > Mx, 
the sfermion masses obtain additional contributions due to SO{10) gauge interactions. These 
increase the values of all cj by the same amount, since all sfermions reside in the 16 of 5*0(10); 
this additional contribution would thus be relatively most important for [31j. However, 
since we need large Higgs representations to realize the breaking chain ([T]), the 5*0(10) gauge 
coupling au hits a Landau pole soon after the unification scale [9]. Hence we expect some 
new, possibly strongly interacting, physics to occur just above Mx- The range of energies 
where 5*0(10) RGE are applicable is therefore probably quite small. 

From Eqs. fl25|) - fl28|l we can derive lower bounds on the ratios of sfermion to gaugino masses. 
Of particular interest for the calculation of the Dark Matter relic density is the relation 

IMiI(Msusy) ^ ^ 

In mSUGRA, the lower bound, which is saturated for ^ "^O' instead slightly below 
unity. This is important, since it implies that for fixed mo and increasing M1/2, one will even- 
tually reach a gaugino mass such that mg^ = fn^o, leading to strong Xi ~ co-annihilation. 
The bound fl29l) implies that this never happens in our scenario. However, as in mSUGRA the 
lighter f mass eigenstate fi can be significantly lighter than e/jlfl We will see later that x? ~ ^1 
co-annihilation remains possible in our model. However, Eq.( l29|) already indicates that the 
parameter space where this can happen is (even) more limited than in mSUGRA. 
Our model also predicts 

^iS^^^^^) > 1 (30) 
|M2(Msusy)| ~ ' ^ ^ 

which means that X2 ^i^d xf decays into SU{2)l doublet sleptons will be strongly suppressed. 
In contrast, in mSUGRA SU{2)l doublet sleptons can be some 15% lighter than SU{2)l 
gauginos. On the other hand, the bound 
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is very similar to that in mSUGRA. It still leaves room for two-body decays of gluinos into 
first or second generation squarks. 



4.3 Mass ratios 

In this section we show numerical results for some (ratios of) masses that are important for the 
determination of the thermal Xi I'^lic density. We focus on masses whose weak-scale values 
are affected by the potentially large Yukawa couplings in the theory. We saw in Sec. 2.1 that 
all third generation Yukawa couplings involving Higgs doublets are quite large at energies 

■l-For this reason, usually the most important co-annihilation channel is Xi—ti co-annihilation [35] : scenarios 
giving TTij^o = mg^ are already excluded, since here fi would be the LSP. 



(a) mg = 1500GeV, = 900GeV, Ag = 0, tanb = 40 (b) = 1500GeV, = 900GeV, Ag = 0, tanb = 40 
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Figure 3: Squared weak-scale running masses in GeV^ of (a) t and f sfermions and (b) the 
up-type Higgs boson, as function of the neutrino mass. The other input parameters are: 
mo = 1.5 TeV, M1/2 = 0.9 TeV, Aq = 0, tan/3 = 40 and > 0. 

> M2] at energies > Mr this includes the new neutrino Yukawa coupling Y^y, which is equal 
to that of the charged lepton by SU{2)ji invariance. The coupling Yj^, which determines the 
Majorana masses of the heavy neutrinos, can also be sizable. Recall that Y^ is related to the 
light neutrino mass and M/j = (a) through Eq. ffT^ . 

Yukawa couplings tend to reduce weak-scale scalar masses for fixed mo and Mi/2- Y^ 
begins to act - on the mass of - at scale M2; at the same scale, the bottom and tau 
couplings become large even if tan/5 is not large, see Eq. (fT2|) . At energies above the 
neutrino coupling Yi, becomes active, reducing the weak-scale masses of fi and of the Higgs 
boson Hu- Above Mc, all weak-scale third generation sfermion masses will be reduced by 
Yn- We therefore expect the difference between first and third generation weak-scale sfermion 
masses to be larger than in mSUGRA. This effect should be strongest for fn and fi- The 
reduction should be more pronounced at small and moderate tan/3, since for large tan/5 all 
third generation Yukawa couplings are sizable even in the MSSM. 

This is illustrated in the left frame of Fig. [3l which shows the dependence of the soft- 
breaking masses of i^, iji, and fji as a function of the mass rrii, of the heaviest light neutrino. 
Recall that this mass is proportional to I/Iat, i.e. smaller rrii, correspond to larger Y^, and 
hence to smaller weak-scale sfermion masses. 

On the other hand, the right frame in Fig. [3] shows that the running soft breaking mass of 
the Higgs bosons with positive hypercharge increases with decreasing m,^. We just saw that 
larger values of Yat reduce m? and m? at all energies below Mx- This reduces the term 
oc Y^. Since this term drives m'^^ to smaller or even negative values, reducing its size leads 
to an increase of the weak-scale value of mj^^. Recall that m|^^(MsusY) < is required to 
achieve electroweak symmetry breaking with tan/? ^ 1. This figure therefore implies that 
the parameter space permitting radiative breaking of the SU{2) x f/(l)y symmetry will be 
smaller for smaller values of m^. 



(a) M,2 = 1200GeV, Ag = 0, tanb = 40, m,=0.4eV 
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(a) = 1200GeV, Aq = 0, tanb = 40, m,=0.2eV 
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Figure 4: The ratio of third generation sfermion masses to the mass of the hghtest neutrahno 
as function of mo, for m,^ = 0.4 eV (left), m,^ = 0.2 eV (right). The values of the other input 
parameters are M1/2 = 1.2 TeV, Aq = 0, tan/3 = 40 and fj, > 0. 

The Yukawa coupling Y^r also reduces the value of the soft breaking mass of S, whose vev is 
responsible for the masses of the heavy neutrinos, and contributes to the SU{2)r x U{\)b-l 
U{1)y breaking. Since this occurs at a scale Mji ^ M^vsy, we introduced a field E, which 
permits to keep the SU{2)r and U{1)b-l -D— terms much below M|.. However, since Yat 7^ 
implies m|(MR) < mlfM^), (S) < (S). Since (E)^ - (E)^ oc (m| - m|) oc M|usy P], 
this effect does not spoil the hierarchy Mr 3> Msusy, but it does give new non- vanishing 
contributions to the masses of sfermions and Higgs bosons. However, they are subdominant 
for most of the parameter space, partly due to the small splitting between Mr and M^, 
and partly because m| and m| receive identical, large gauge contributions, in particular for 
Q > Mq where S and S are embedded in (anti-)decuplets of SU{4:)c- In fact, we found that 
the ratio (m| — m|)/m| is, at most, a few %. Hence, these new D— term contributions can 
be ignored. 

Fig. m shows the dependence of the ratios '^i^ rI'^xi ^^"^ bI'^xIi taken at scale Q = 
^susY, on niQ. For ml -C M^^^ the stop squarks are significantly heavier than the stau 
sleptons. This qualitative behavior is the same as for first and second generation squarks, see 
Eq. fl28|) . On the other hand, if mo > M1/2 and relatively small (left frame), m^^ can be 
smaller than m^^^, since the top Yukawa coupling is significantly larger than that of the r 
lepton. 

The right frame of Fig. H] shows that increasing Yj^ reduces the dependence of third gen- 
eration sfermion masses on mo. In fact, for m^ = 0.2 eV we observe a sort of "focus point" 
[37] for mi^, i.e. m^^ (Msusy) becomes almost independent of mo. This implies that there 
is no focus point behavior of m^^ (Msusy), i-e. this soft breaking parameter, which largely 
determines electroweak symmetry breaking for tan^/? ^ 1, does depend on mo. Hence large 
values of mo will not be "natural" (by the definition employed in refs.|37]) if Y^ affects the 
weak-scale third generation masses significantly. 

The scalar masses shown in Fig. [Hare running masses at scale Msusy- The physical masses 




Figure 5: The ratio of m^-;^ to m^o (left), m^o to 2m^o (right), as the unification scale Mx is 
varied. Recall that this implies corresponding variations of the intermediate scales Mq and 
Mr, see Fig. [D 



will be affected by threshold corrections and, more importantly for third generation sfermions, 
by mixing between SU{2) singlets and doublets. This mixing reduces the mass of the lighter 
eigenstates fi and ti, so that ruf^ < min(mf^, rrifj^) and similar for m^^. Nevertheless Fig. H] 
shows that co-annihilation will usually only be possible with fi. This is similar to mSUGRA. 
Note, however, that Fig. |4]is for tan/5 = 40. Recall that co-annihilation with first or second 
generation sfermions is not possible here, see Eg. (1291) . Moreover, comparison of the two frames 
of Fig. m shows that the effects of Yn on the f masses are quite small if mg <^ ^1/2- -^^ ^ 
result, we find that f co-annihilation is possible in our model only for tan/5 > 27. 

The left frame of FigjS] illustrates the dependence of the ratio mf^/m-j^o on the GUT scale. 
Recall that for Mx = 10^^'^ GeV our model becomes indistinguishable from mSUGRA, as 
far as the weak-scale spectrum is concerned. As discussed in the previous Subsection, in the 
absence of new large Yukawa couplings this ratio can only become larger as the intermediate 
scale is turned on. However, we saw in Fig H] that the Majorana Yukawa coupling Yx does give 
a large positive contribution to the RGE of rrif^, reducing its weak-scale value. These two 
effects clearly compete with each other. We see that even a rather large Yx, corresponding 
to m^^ = 0.2 eV, can change ^^^(Msusy) significantly only if Mx < 10^^'^ GeV; recall from 
Fig. [1] that this corresponds to the region of parameter space where Mr < Mc- In this case 
the possibility to have fi co-annihilation obviously strongly depends on m^^. 

The right frame of Fig I5] shows the ratio mAo/2mj^o. This ratio needs to be close to unity 
for Xi annihilation through s— channel exchange to be enhanced. We see that reducing 
Mx, i.e. turning on the intermediate scales, slightly increases this ratio even if Yx is small. 
For given mo, this can be compensated by increasing Mi/2- We thus expect the "A— funnel" 
region to survive in our scenario, if tan/5 > 50 and for small Yx- Increasing Yx will increase 
myi(MsusY); this is analogous to the increase of depicted in the right frame of Fig. [31 




Figure 6: Constraints on the (mo, M1/2) plane of our model. The grey areas are those excluded 
by the EWSB condition or by tachyonic or too light sfermions. The region excluded by the 
Higgs and chargino mass constraints is shown in bright red and the own excluded by the 
b s'-f constraint in pink. The blue area satisfies the — 2 constraint fl20l) . while green 
regions satisfy the Dark Matter constraint ( !2T|) . Finally, black regions satisfy all constraints. 



4.4 Regions of the (mo, M1/2) plane 

In this Subsection, we show the (mo, M1/2) plane of our model, indicating the regions where 
the various accelerator as well cosmological constraints discussed in Sec. 3 are satisfied. We 
scan the parameter space only up to (mo, M1/2) = (2000 GeV, 1500 GeV). Even larger sparticle 
masses appear quite unnatural. The LHC should be able to probe the entire parameter space 
we show [38]; recall that M1/2 = 1.5 TeV corresponds to a gluino mass around 2 TeV in our 
scenario. We focus on large values of tan/?. We saw in the previous subsection that this is 
required both for the A— funnel and the for f co-annihilation region in our scenario. Finally, 
sign(/i) is chosen positive in all plots, in accordance with the indication of an additional 
positive contribution to (7^; recall also that taking /i > makes it easier to satisfy the 6 — 57 
constraint [27] . 

A first example, for Aq = and tan (3 = 40, is presented in Fig. (6] the left (right) frame 
is for small (large) coupling Y^. The grey regions are mostly excluded by the requirement of 
correct electroweak symmetry breaking; in the right frame the region of small mo is instead 
excluded because fi is too light (below either the LEP limit or the mass of x?)- As expected 
from the discussion of Fig. [3l this region is considerably larger for large Y^- 

The bright red regions are excluded by the chargino search limit (fT6!) or by the limit (fT5|) 
on the mass of the lightest CP-even Higgs boson; the latter is relevant for M1/2 ^ 500 GeV, 
while the former excludes the narrow red strip bordering the grey region at large mo and large 
Ml/2- Finally, the pink regions are excluded by the constraint fll7p on the branching ratio for 
radiative b decays. Some supersymmetric contributions to the corresponding amplitude grow 
oc tan/?. This constraint therefore becomes relevant at the large values of tan/? required to 
realize f co-annihilation and/or the A— funnel in our model. 



Turning to observables that require a non-vanishing contribution from supersymmetric 
particles, in the blue regions the constraint fl20l) from the anomalous magnetic moment of 
the muon is satisfied. The corresponding diagrams are quite similar to those contributing to 
6 — *• 57 decays. In particular, some contributions again grow oc tan/5. As in mSUGRA [14j, 
we find regions of the parameter space at sufficiently large M1/2 where electroweak gauginos 
and sleptons are sufficiently light to give a sizable positive contribution to g^, while (stop) 
squarks are sufficiently heavy not to reduce the branching ratio for b s'y decays too much. 

Note that the red, pink and blue regions all extend to much larger values of M1/2 than in 
mSUGRA [13]. The reason is that the corresponding constraints probe weak-scale (s)particle 
masses; we saw in Eqs.( 12^ - (I25]) that a given M1/2 corresponds to much lighter gauginos and 
sfermions in our scenario than in mSUGRA. Moreover, we saw in Figs. [3] and H] that the 
additional large Yukawa couplings in our model tend to reduce weak-scale stop masses. They 
also increase fn'jj^, which leads to a reduction of |yu| via the condition of electroweak symmetry 
breaking. Both effects, which become more important for smaller m,^, increase the absolute 
size of the stop-chargino loop contribution to b —>■ s'j decays. This has to be compensated by 
increasing ttiq and/or Mi/2- The b sj constraint is therefore relatively more important in 
our scenario than in mSUGRA, especially if is sizable. 

Note also that the region excluded because it does not permit radiative symmetry breaking 
has a pronounced slope even for the larger neutrino mass, i.e. smaller coupling Y^. This 
shows that Tn'j^^[MsvsY) has significant dependence on mo, as remarked earlier, i.e. there is 
no focusing behavior of this parameter. As expected from our discussion of Fig. [3l this upper 
bound on mo becomes stronger when Y^ is increased, i.e. when m^, is decreased. In a strip 
close to this excluded region we nevertheless expect the lightest neutralino to have a large, 
perhaps dominant, higgsino component; this region will therefore have a somewhat similar 
phenomenology as the "focus point" region in mSUGRA [37j, especially as far as Dark Matter 
is concerned. 

Finally, in the narrow green strips the constraint (pTj) on the Dark Matter relic density is 
satisfied; these strips would obviously look broader if we had indicated the 2a allowed region, 
as more commonly done. The overlap between the DM- and (7^— allowed regions is colored in 
black. 

In Fig. |6] we find two such regions. At small mo Xi is bino-like, and achieves a sufficiently 
small relic density through co-annihilation with fi. For small Yat (left frame) this region is 
strongly constrained by the bound on b —>■ sj decays. We saw in Fig. H] that increasing Y/v 
reduces the f masses, making it possible to find scenarios with rrif^ ~ m^o even if M1/2 is 
large. 

We just saw that for values of mo not far below the upper bound imposed by electroweak 
symmetry breaking, Xi has a sizable higgsino component. For some range of parameters it 
achieves the correct relic density mostly through annihilation into channels involving weak 
gauge bosons. As in mSUGRA, this second DM-allowed region extends to very large values 
of mo and M1/2, with x^ becoming increasingly higgsino-like (and therefore co-annihilation 
with X2 ^^"i xf becoming increasingly important [39j.) 

As in mSUGRA, tan (3 = 40 is not large enough to allow rriA — 2m^o if /i > 0. Fig. [7] 
shows that this "A— pole" region becomes accessible for tan (3 = 50. Sufficiently small values 
of niA are only possible if the soft breaking mass mfj of the second Higgs boson also becomes 
negative (and large) at the weak scale. We saw in the discussion of Fig. [5] that decreasing 
rriu will increase m|^ (Msusy)- Indeed, in Fig. [7] we find a well-defined A— funnel only for 
nil, = 0.4 eV (left frame). 




Figure 7: Constraints on the (mo, M1/2) plane of our model. Parameter values and color code 
are the same as in Fig. [6l except that tan (3 has been increased to 50. 
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Figure 8: Number of allowed points (by all constraints) for ITeV < mo < l.STeV; l.lTeV 
< Ml/2 < 1.4TeV, with the grid 25GeV. The allowed region is very large when tan/3 = 49. 



If we instead take rrti, = 0.2 eV (right frame), we find that the Xi relic density becomes 
too low in the entire allowed region of the (mo, M1/2) plane we scanned. One reason is that 
increasing reduces /^(Msusy), as discussed above. This increases the coupling of the LSP 
to neutral Higgs bosons, in particular to A. Since for tan /3 = 50 the b and r Yukawa couplings 
are quite sizable, virtual A exchange diagrams become large, even though 2m^o is somewhat 
below rriA- Increasing mo increases tua, but at the same time decreases /i even further, and 



therefore does not allow to achieve a DM relic density above the lower bound in the range 
(12T1) . Moreover, recall that reducing m^, also reduces the f masses. In addition, the very large 
value of tan (3 considered in this figure leads to large fifji mixing, which allows XiXi ~^ t^t^ 
annihilation through f exchange even if the initial state is in an 5— wave |40j. Finally, for 
Ml/2 close to its lower bound, x^fi co-annihilation again becomes important. Note that this 
indicates that the DM-allowed region may be quite large for some tan/3 between 40 and 50, 
and rriy = 0.2 eV. Indeed, Fig. [8] shows that for tan/5 = 49, about 50% of the points we 
scanned that satisfy the other constrains are also compatible with the DM constraint. 



(a) tanb=40, Ao=mo, my=0.4eV (a) tanb=40, Ao=mg, m^=0.2eV 
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Figure 9: Constraints on the (mo, M1/2) plane of our model. Parameter values and color code 
are the same as in Fig. El except that we now take Aq = uiq {Aq = —mo) in the top (bottom) 
row. 

In Fig. iniwe explore the effect of taking a non-zero value of Aq. We see that the value of Aq 
can have quite a dramatic effect on the region excluded because it does not allow electroweak 



symmetry breaking. This can be understood as follows. By dimensional arguments and the 
fact that scalar masses always appear as squares in the RGE, the soft breaking mass of the 
up-type Higgs boson at the weak scale can be written as 

mjj^ (MsusY ) = aml + hMl/^ + cAl + dM^/^Ao . (32) 

The values of the coefficients a, 6, c, d depend on the dimensionless couplings in the theory, 
as well as (logarithmically) on Msusy- In our model, a and d are positive while h and c are 
negative. Hence increasing mo makes EWSB more difficult, while increasing M1/2 makes it 
easier if M1/2 > \Aq\. This explains the qualitative feature of the regions excluded by the 
EWSB constraint in Figs. [6] and [71 

On the other hand, if |Ao| ^ M1/2, increasing the absolute value of Aq also aids EWSB 
independent of its sign. This explains why the region excluded by the EWSB constraint 
becomes much smaller in the two left frames of Fig. O Fig. E] shows that, for the given small 
value of Yat, the EWSB constraint only excludes scenarios with uiq > M1/2 even if = 0. 
In the critical region |Ao| = ^0 is thus always sufficiently larger than Mi/2- Finally, for given 
absolute value of Aq, EWSB will be easier for negative than for positive Ao. This explains 
why the EWSB excluded region is significantly larger in the upper-right frame of Fig. [9] than 
in the lower-right frame. Note also that a sizable decreases the absolute size of c, since 
Yn reduces \At\ for Q > Mc, see Eq.(A.34) 

A nonvanishing Aq also changes the regions allowed by the other constraints. In particular, 
Aq < increase ti — tn mixing. This has two effects. On the one hand, it increases the 
radiatively corrected mass of the lightest CP-even Higgs boson, thereby reducing the size of 
the red regions in Fig. O On the other hand, it increases the tx^ contributions to radiative 
b ^ s'y decays, increasing the size of the pink regions. This latter effect completely removes 
the DM-allowed region close to the EWSB-forbidden region, where Xi h^is sizable higgsino 
component. As a result, for Aq = —mo, only the small f co-annihilation region survives. 
On the other hand, for Aq = mo we again find sizable DM-allowed regions at large mo; the 
structure in this (black) region at M1/2 — 800 GeV in the top-left frame is due to the opening 
of the XiXi ~^ tt channel. 

5 Summary and Conclusions 

Supersymmetric 5*0(10) GUTs have become attractive extensions of the SM, especially since 
the observation of the nonzero neutrino mass. However, there is a small discrepancy between 
the order of the expected mass of the right handed neutrino in the seesaw mechanism and the 
GUT scale; it can be explained in a natural way if one postulates intermediate scales where 
the gauge symmetry is larger than that of the SM, but smaller than 5*0 (10). 

Therefore, in this work we chose a model [0] which gives us intermediate symmetry breaking 
scale(s), and analyzed how this affects the low energy phenomenology. We found that the 
relation between weak-scale and GUT-scale parameters is quite different in this model than 
in the widely considered mSUGRA scenario. Perhaps more importantly, ratios of different 
weak-scale masses also differ from mSUGRA. In particular, the slepton to electroweak gaugino 
mass ratios are higher than in mSUGRA. As a result, co-annihilation is only possible with 
the lighter f eigenstate, and only at large tan/? and/or large Yukawa coupling Y/v of the SM 
singlet neutrinos; the latter corresponds to small values for the light neutrino masses. 



Radiative electroweak symmetry breaking also is more difficult in this model than in 
mSUGRA. This makes it easier to find Dark Matter allowed solutions where the lightest 
neutralino has a significant higgsino component. As in mSUGRA, the location of this region 
strongly depends on Aq; in addition, we find a strong dependence on Y]y, i.e. on the light 
neutrino mass. We also found that for very large tan/3 and large Yj^ most of the (mo,Mi/2) 
plane leads to too small a Xi relic density. As a corollary, there exist combinations of Yn and 
tan/? where ^dm^^ has weak dependence on mo and M1/2; however, in this case it depends 
strongly on tan/5 and Y^. Finally, as in mSUGRA the A— pole region only exists at large 
tan/3; it disappears for large values of Yat. 

We would like to point it out that, even though our analysis is done for a specific model, 
many of our results should remain qualitatively correct for other 5*0(10) GUT scenarios, as 
long as the seesaw mechanism at an intermediate scale plays a role. In particular, the relation 
between the right-handed stau mass and the Majorana Yukawa coupling Y^, which largely 
determines the behavior of the co-annihilation region, does not depend on the details of either 
the symmetry breaking chain or the seesaw structure. Any partial unification above the see- 
saw scale also implies that Y/v will affect other sfermion masses, and hence the conditions for 
radiative symmetry breaking. 

In summary, the model we considered relates several phenomena, and can hence be probed 
through a large variety of measurements, from proton decay (which imposes limits on the 
GUT scale) over neutrino masses and Dark Matter physics to collider physics. We intend to 
investigate characteristic features of this scenario at the LHC in a future publication. 
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A Renormalization Group Equations 

In this section we list all relevant one-loop renormalization group equations explicitly. Our 
calculations are based on the general expressions of ref . . We divide the entire energy range 
between the SUSY and GUT scales into five regions, with different particles participating in 
the RGE and different symmetry groups: 

• Region I (MsusY < Q < M2) : SU{3)c x SU{2)l x [/(l)y 

• Region II {M2 < Q < Mn) ■ 5f/(3)c x SU{2)l x f/(l)y 

• Region III {Mr < Q < M^) : SU{3)c x U{1)b-l x SU{2)l x SU{2)r 

• Region IV {M^ < Q < Mc) : 5f/(3)c x U{1)b-l x SU{2)l x SU{2)r 

• Region Y {Mc < Q < Mx) ■ SU{A)c x SU{2)l x SU{2)r x D 

In the following Subsections we discuss the running of the supersymmetric parameters 
(gauge couplings and parameters of the superpotential) and of the soft breaking parameters, 
respectively. 



A.l Superpotential Parameters 



We begin with the parameters that preserve supersymmetry. The running of the gauge cou- 
phngs is described by 



l-^Pg^gl with = S{R) - 3Ca{G) . (A.l) 



d 

R 



Here t = hi{Q/Qo), a labels the factor group, R the representation of the matter and Higgs 
superfields under this group, Ca is the quadratic Casimir of this group, and the Dynkin index 
S{R) is defined by Tr(t^t5) = S{R)6ab, tA,B being matrix representations of the gauge group. 
Our notation for a generic superpotential is 

w = h^'^'^^i^.^u + ■ (A.2) 



The running of the parameters appearing in ( 1A.2I) is given by 
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li + (j ^ i) , (A.3) 

where summation over repeated indices is understood. The anomalous dimensions 7^ are given 

by 

li = \Y,,,Y^^^-25lglC^{i). (A.4) 

In our case, the superpotential below Mx has been given in Eq.Q for region I, in Eq.Q 
for region II, in Eq.(l8]) for regions III and IV, and in Eq.(l9]) for region V. Recall that we take 
Y2 = Yq^2 = Yi^2 = Yu^2 = Yfifl = Ye^2 = 0; for the sake of simplicity we therefore suppress the 
superscript 1 on the Yukawa couplings in the following. These couplings are 3x3 matrices 
in generation space. We will write the RGE for general matrices, although we only kept third 
generation couplings in our numerical analysis. We use the general notation 

Jt^ = 1^^- (^-^^ 

where / stands for any matter fermion. In the following we list these as well as the gauge 
beta-functions in the five different energy regions. 

A. 1.1 Region I 

The coefficients of the gauge beta functions are 

(3g^ = (33/5, 1, -3) for a = (ly, 2^, 3c) , (A.6) 



where we have used GUT normahzation for the U{1)y factor. The corresponding coefficients 
for the MSSM Yukawa couphngs are 

Py. = >;(7l + 7g) + ; 



where 



7| = 

7f = -^eY^'-To'l-lal, 



7S = tr{3Y,Y^ + Y^Yj)-^gl-lgl; 

= ^trY^Y^-^gf-lgl (A.8) 



A. 1.2 Region II 



= (12, 2, -3) for a = (ly, 2l, 3c) . (A.9) 

The Yukawa couphng beta functions of the MSSM matter fields have the same form as in 
Region I, but we need to introduce an RGE for Y^: 

Py^ = 1^7v(7l + 7I) + iEYn . (A.IO) 

Except for 7^ the anomalous dimensions of the MSSM matter fields also remain form invariant, 
and we have to introduce an anomalous dimension for 6 : 

7I = 2y;ye + 4>V - ^^7? ; 
1 94 

7I = ^tT{Y},Y,,) - -gl (A.11) 

Recall that we are now dealing with the couplings !}■ 1 (/ = u, d, e), which are related to the 
MSSM couplings Yf via Eqs.([l0])-([l2]). 

A. 1.3 Region III 



Pg^ = (15, 2, 6, -3) for a = {1b-l, 2l, 2r, 3c) , 



(A.12) 



where we have again used GUT normahzation for the U{1) couphng. The effective coefficient 

^ -6 + ^ 

5 ^ ^ 5 



48/5 for the running U{1)y couphng hsted in Table 2is|-6 + |-15, which foUows from the 



matching condition gy'^ = ^g^^ + 19b-l- 

Since the underlying symmetry group is enhanced, and the matter superfields form multi- 
plets correspondingly, their anomalous dimensions receive contributions from the heavy gauge 
bosons that become active in this energy range. We switch to the notation of Eq.Q, i.e. we 
introduce Q'^ instead of If^ and D'^, while E'^ and A^^ are united in L'^ and Hu and are 
united in The number of independent Yukawa couplings is thus reduced to three: 

Py, = n(7g: + 7l)+7gn; 

Py^ = YNh^ + ll)+l^YM. (A.13) 
The relevant anomalous dimensions read: 



7f: 


= 2Y^Yi + h},Y^-'^-gl- 


3 2 

'^9b-l 


7f 


= 2Y,Y^-\gl-\gl^,- 






= 2YqY^ — -g2 — -^Qb-l ~ 


8 2 

- 3^?3. 


Q" 
7q^ 


= 2Y^Yq — -g\ — —gs-L 


8 2 


it 


= tr(3F,i;t + YiF/) - ^(7^ 


3 2 

2^72, 


7| 


= hriYj,Yi,)-Agl-3gl 


-L ■ 



(A. 14) 

Here we have continued to use g2 for the SU{2)l coupling, and denoted the SU{2)fi coupling 
with gii. Eqs.f lA.Mj) are consistent with [42j, taking the appropriate normalization. 



A. 1.4 Region IV 



= (15, 2, 6, 0) fora = (lB-L,2L,2fi,3c). (A.15) 

Since the new massive fields becoming active in this energy range are singlets under U{1)b-l x 
SU{2)r X SU{2)l, only the running of the SU{3)c group changes. Moreover, the Yukawa 
coupling beta functions are those of Region III. 



A. 1.5 Region V 



= (42, 42, 34) for a = {2l, 2r, 4c) • (A.16) 

Since many new fields become active at Q > Mc, all gauge /3— functions increase quite dra- 
matically. In GUT normalization, gs-L = Qi = 94,, where is the S'f/(4)c gauge coupling; 
this explains the entries in the last row of Table 2, with = | ■ 42 + | ■ 34. 

No new Yukawa couplings appear in this energy range; instead, the couplings Yg and Yi 
get unified into the single coupling Y. At the same time, all MSSM matter superfields are 



now in F or F'^ introduced in Eq.Q; the D symmetry ensures that the anomalous dimensions 
of these two superfields are the same. Moreover, the SU{2)ji triplet Higgs superfield 5 gets 
embedded into the much larger representation S/j, and the partner also appears, with 
identical anomalous dimensions. 
Therefore Eqs. flA.13p change to 

Py. = y^(7l + 7g)+7l>^^. (A.17) 
The anomalous dimensions appearing in Eqs. ( lA. 171) are: 

7| = tr(4rrt) _ _ 3^2 . 
7g = hTiY^Yi,)~Agl-9gl (A.18) 

A. 2 Soft SUSY-breaking Parameters 

We write the part of the Lagrangian that softly breaks supersymmetry as 

CsB = - lb''<P^<P, - - ^M^XAa + h.c. (A.19) 

We assume universal boundary conditions, 

{w?)l = ml5l ; 

Ma = My2^a, (A.20) 

which hold at scale Q = Mx- Here Y^^^ are the superpotential couplings introduced in 
Eq.(IS21). 

The /5— functions of the soft breaking parameters are defined by 
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PI 

2^'- 

'^gXPg^ . (A.21) 

Here Pg^ are the coefficients of the gauge /3— functions introduced in Eq. flA.ll) . The other 



/3— functions appearing in Eqs. flA.2ip can be written a£| 



1 



1,, o.. 1, 



85^M.Mt^^^,(^) + 2^,2(t^)^tr(t V) . 



(A.22) 



The last term in /3„2 can be nonzero only for f/(l) group factors. In the case at hand it is 
therefore either proportional to 



O 2 2 . L \ 2 2o2|2| 



(A.23) 



or to 



Sb-l = o (6"^^ — 6m| + tr[2mQ — 2mQc — 2m^ + 2m^c]) . 



(A.24) 



For better readability, in Eqs. (]A.23P and flA.24p . as well as in subsequent equations, we have 
omitted the tildes on the subscripts of the squared scalar soft breaking masses; moreover, we 
use ms and for the soft mass of whatever parts of the original S and S superfields are 
active in a given energy range. Note that both Sy and Sb-l evolve homogeneously. Since the 
boundary condition flA.20p for scalar soft breaking masses implies Sy = Sb-l = at scale 
Mx, they vanish at all scales. For completeness we nevertheless list these contributions in the 
following. 

We are now ready to give explicit expressions for the soft breaking /5— functions in the 
energy regions defined above. 



A. 2.1 Region I 

Here the RGB are those of the MSSM HSl: 



= K 

+ Yu 

= hd 

+ Yd 

= K 

+ n 



16 



13 



tr(3F„yJ) + Sy.Fj + YjYd - -fgl - 3gl - -^gf 

tr(6/i„rj) + Ah^Y^ + 2Y}hd + ^dlMs + Gg^M^ + '^gfM, 

16 7 
tr(3r,ri + FeF;) + 5YdYj + Y^ " -jdl - ^qI - ^dl 

32 14 
tr(6/i,r] + 2KY}) + AYlhd + 2Y^K + ^qIM, + QglM^ + —glM, 

tr(3y,ri + Y,Y}) + hY}Y, - 3^72^ - ^1 

tT{6hdY} + 2KY^) + + QglM2 + y^/^Mi 



(A.25) 



*We suppress terms that can be nonzero only in the presence of complete gauge singlet chiral superfields. 



Pr, 



Pml 
Pml 



Pe 



= B 



tr(3i;Fj + 3F,rJ + Y,Y}) - 2,gl - -gl 

tr(6/i„rj + QhdYi + 2heY}) + &glM2 + \glMi 

o 



(A.26) 



tr 



6 



{ml + 2m|jy„ty„ + {ml + 2m^jyiy, + [yjn + Y}Y,]ml + 2y>^^„ 



652'|M2|^--^i^|Mi| 



2 2 c 



32 



+ 2YlmlYa + 2/it/i, + 2/i>, - -gl\M^\' - Gg^lM^l' - -glm' - -glSy ; 



K + 2mijy;ye + 2y;m^ye + r^ml + 2/it/i, _ 6^2^171^2!^ - -gl\M,\^ - -glSY ; 

32 32 

(2m^ + 4m|Ji;Fj + 4F„m|Fj + 2F„Fjm^ + - -^gl\M,\' - -gl\M,\^ 



:9lSy ; 



32 



(5^. = {2ml-rAm]j^)YaYl-rAYamlYl-r2Y,Ylml-rAhah\--gl\M,\^--^^^^ 



+ 



24 6 

(3ml = (2m^ + 4m|jyen^ + 4yemiy; + 2yey;m^ + 4/ie/iI-^^i|Mi|2 + -^2^y. (A.27) 



A.2.2 Region II 

Most expression from Region I remain form-invariant; however, the Yukawa couphngs should 
now be interpreted as the high-scale couplings Yf i rather than as low-scale (MSSM) couplings 
Yf. In addition, the beta-function for the SU{2)l singlet slepton mass changes, and we have 
to introduce beta-function for m^ as well as Hn- 

1 36 
Ph^ = hN[^tT{Y^YM)+4YeYj + 2YNYi,--gl] 

+ yjv[y>iv + 8yX + 4y>iv + y^?|Mi|2]. (a.28) 

= £^[itr(y^y^)-^^?]^Ms[trO(y>^) + ^5i|Mir]. (A.29) 

(3^2 = {2ml + 4mj,jY,Y^ + AY,mlY^ + 2Y,Y^ml + mlYNYl + 2YlmlYN 
+ Y^Y^ml + 2Y}^m^YN + (4m^ + 2m|)y7vy^ + 4/ie/ii + 2/ijv/i)v 

tr[^y^y^m| + 2YlmlYN + lY^Y^ml] + h^h^^ - ^gl\M,\' + ^glSY{A.30) 



A 



-I 



Ms appearing in the Eq. flA.29P is the supersymmetric S and S mass, which comes from a 
term M^SS in the superpotentiaL Note that to one-loop order By. does not appear on the 
right-hand side of any other RGE, hence it has no impact on the low-energy spectrum. We 
nevertheless list its RGE for completeness; it might be relevant, e.g., for the detailed dynamics 
of intermediate-scale symmetry breaking, which we here merely parameterize through the vev 
a. 



A. 2. 3 Regions III and IV 



Here the number of independent parameters diminishes: S'f/(2)ij invariance implies = 
rriu = mQc, rrie = rri]^ = ttilc, mu^ = rriH^ = m$, hu = = hg and he = hj\f = hi at energies 
>Mr: 



+ y, 

= hi 

+ Yi 

= hN 



+ Yi 



N 



32 

tr(6/i,Fj + 2hiYi^) + AhgYg^ + 2Yg^h, + —glM^ + QglM^ + QglM^ + 
tr(3F,F; + YiY,^) + Sl^F/ + Y,Wi + h^l^Y^ - ?,gl - 3^?| - ^-gl^^ 
tr(6/i,r; + 2hiY^) + AhiY^ + 2Y^hi + 3F>;v + Q9IM2 + QglM^ + 2g 

itr(r^r^) + AYiY,^ + 3r^r^ - 7gl - \gl^^ 

y>iv + ^Y^hi + 6F>^ + UglMR + 9^|_^ 



M 



B-L 



B-L 



(A.31) 



Pb 



B 

By 



tr(6r,F; + 2YiY^) - ?,gl - 3g 
MYnY^) - 8gl - 6gl_^ 



+ /i 

My 



tr(12/i,F; + 4/i,F/) + QglM2 + GglMn 
tr(F>^) + WglMn + 12(7|_^Mb_l 



(A.32) 



(3^2 = tr [6(m| + rjr, + GY^ml^Yg + 2(m| + mi)r/r, + 2r/mier/ 



^-1 



P. 



+6hlhg + 2h]h, 



6gl\M,\' - 6gl\Mn\' ; 



16|Mk|2^2 _ i2\Mb-l\'91_l + S^Ll-^b-l ; 
+ 4hlh,--gl\Ms\'-6glm' 

= 2(mi + 2ml.)r/r, + 2r/r,mi + 4r/m^y, 

+ 4/i;/i,-6^2|^2r-3^|_jM, 

= (2m^. + 4m|)F,r; + AY^mlY^ + 2F,F>Je + 4/i,/i; 

1 



ib-l\ — -^9b-l^b-l ; 



(A.33) 



32 1 



(2mi^ + 4m|)F;F/ + 4F;miF/ + 2YiYi^ml. + Ahihj + -ml^Y^Y^ + -Fiv^mi. 



+ 3yjimicyiv + 3yjim|yjv + 3/ijv/iiv - Q91\Mr\^ - ^gl_L\MB-L? + ■ 

Note that E, which we introduced to allow a D— fiat direction for symmetry breaking at scale 
Mr, does not have any superpotential couplings, hence its soft mass does not appear in any 
of the other RGE. We again hst its RGE for completeness. 

A.2.4 Region V 

At scales above Mc the spectrum further simplifies: 6*422 invariance implies that mq — uil = 
mp-i rriQc — rriLc — m^c, and hg — hi = h. In addition, the discrete D symmetry implies 
m^^ — = and = m^c: 



Ph ^ h 
+ Y 

+ Y^ 



tr(8/iFl") + —HnYI + 4hY^ + 2^1" /i + Ibg'lM^ + 6c/^M2 + G^lM^ 



^tr(F^F^) + 4yyt + l^Y^Y^ - ^gl - 7g 
trihNY},) + 8Y^h + 15y>^ + 335^M4 + Ug^Mn 



(A.34) 



Pb ^ B [tr(8FFt) _ 3^2 _ 3^^] ^ ^ [tr(16/irt) + Qg^M^ + GglMn] ; 



-tr(y;vli ) - - 18^4' 



+ ms 



tr(y>jv) + 16glMR + 36glMA .(A.35) 



13ml 


= tr [8(m| + ml)Y^Y + SY^mlY + Shlhu] - QglW^l^ - 6(/||Mb|2 ; 
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= Ui{YlYMm^ + r^r^m|) + tr(24m|F;v + /^U^) - 16|M^|2^| - 36|M4| 






= -lQ\MR\^gl-?,Q\M,\^gl- 






1 ^ 1 c; 
= 2{ml + 2m|)Fty + 2Fty^^ + + 2m|)Fjv>l + — F^F^m| 






15 15 

+ — + —hNh\ + 4FWFF + 4.h^h - 15(7^|M4p - 6^72^2!' ■ 


(A.36) 
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